Deformed Fermi surfaces in ultracold Fermi gases 
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The superfluid transition in an ultracold two-component atomic Fermi gas is analyzed in the 
case where the two components have different densities. We describe a superfluid state which 
spontaneously breaks the rotational-symmetry by deforming the Fermi surfaces of both species into 
ellipsoidal form. At relatively large hyperfine-spin asymmetries, this deformation is shown to help 
the appearance of pairing, which in the rotationally-symmetric (BCS) case would be forbidden by 
Pauli blocking. The prospects for experimental detection of such a deformed Fermi surface phase 
are discussed. 



S3 ■ 
O ■ 

o ; 

u • 

Or 

3 ; 



-a 
c 
o 

(N 
> 

in 
o 
o 

-I— > 



o 
o 



X 



The present capabilities of cooling of atomic ensem- 
bles allow for reasonable expectations to observe a su- 
perfluid transition in ultracold fermionic systems, in di- 
rect analogy to the BCS superconductivity [lj . The tem- 
peratures achieved in recent experiments on fermionic 
atoms are a fraction (~ 0.1 — 0.3) of the Fermi- 
temperature, i.e. atoms in a trap are in the quantum 
degenerate regime and, therefore, attractive two-body 
forces are expected to drive the Cooper instability. The 
strength of the two-body interactions can be tuned us- 
ing a Feshbach resonance by varying the external mag- 
netic field thus the entire range from weak to strong 
couplings can be probed. In the crossover region the 
Feshbach resonance may strongly enhance the pairing 
interaction and give rise to high temperature superflu- 
idity 0. Recent experiments have probed the conden- 
sation of fermionic pairs above the Feshbach resonance, 
where the system does not support a genuine two-body 
bound state H, IE ■ Whether the observed correlated 
pairs are, in fact, weakly bound and spatially extended 
Cooper pairs, is not clear yet; however the measured col- 
lective modes of 6 Li atoms under these conditions are 
consistent with the superfluid hydrodynamics of a Fermi- 
gas and provide evidence for superfluidity in a resonantly 
interacting Fermi gas 0. 

The very low temperatures (in the nanokelvin range) 
and densities reached in the experiments considerably re- 
duce the contribution from L ^ orbital angular mo- 
mentum waves to atomic collisions. Therefore, s-wave 
collisions, which can be characterized by the scattering 
length a, are the most relevant for the description of these 
systems. As usual, we take a < to indicate an attrac- 
tive interaction between the atoms. Since Pauli's exclu- 
sion principle forbids s-wave interaction between indis- 
tinguishable fermions the pairing should appear between 
fermionic atoms belonging to different hyperfine states. 



Such systems, where two hyperfine levels arc populated, 
have been created and studied experimentally with 6 Li 
and 40 K atoms 0, E 0. The BCS theory predicts a 
suppression of the pairing correlations when the Fermi- 
energies, or equivalently, the densities of the two hyper- 
fine states |1) and |2) are different. In the low density 
limit kp\a\ <C 1 (kp is the Fermi-momentum) the value 
of the critical density asymmetry a — (pi — p<x)j(j)\ + p%), 
where pi } 2 are the densities of hyperfine states |1) and |2), 
for which the superfluidity vanishes can be deduced ana- 
lytically |8| • The dependence of the pairing gap A at the 
Fermi surface on the total density p = pi+p2 = kp/ (3n 2 ) 
and the density asymmetry a is described by 
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where Ao ~ 8e _2 /iexp [— 7r/(2fcir|a|)] C 1 is the gap 
in the symmetric matter and p is the chemical po- 
tential. Therefore the gap disappears for asymmetries 
a > aimax = 3Ao/(4/i), which in this limit is a very 
small number. For example, for the pairing of 6 Li 
atoms in the states |1) = \F = 3/2, nip = 3/2) and 
\2) = 3/2, 1/2) , for which the triplet scattering length 
is a = — 2160as (where as is the Bohr radius), and 
at the density p — 3.8 x 10 12 cm~ 3 (corresponding to 
kp\a\ = 0.55) the maximum asymmetry at which BCS 
pairing is possible is only a max = 0.07. 

The purpose of this work is to demonstrate that the 
superfluid state in ultracold atomic gases can persist for 
density asymmetries a > a max , and can be enhanced in a 
range of a < a max due to spontaneous deformation of the 
Fermi spheres of two-hyperfine states in the momentum 
space. It has been shown earlier (in other contexts) that 
the deformation of the Fermi surfaces of asymmetric two- 
component superconductors into ellipsoidal form leads 
to a novel ground state with deformed Fermi surfaces; 
the associated superconducting phase brakes global rota- 
tional symmetry of the space from 0(3) down to 0(2) 0. 
The deformed Fermi surface superfluid (DFS) phase be- 
longs to the class of the superconducting states with bro- 
ken space symmetries, an example of which is the Larkin- 
Ovchinnikov-Fulde-Ferrell (LOFF) phase 0; which has 
been studied in ultracold atomic gases in Ref. 01 ■ The 
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underlying principle that makes these phases favorable is 
the compensation in the increase in the kinetic energy of 
the superconducting phase by the gain in the (negative) 
condensation energy due to a rearrangement of quasipar- 
ticle distributions. In the LOFF phase the compensation 
is achieved by sampling Cooper pairs with finite total 
center-of-mass momentum; in the DFS phase the same 
is achieved by a deformation of Fermi surfaces at zero 
total momentum of the pairs. Formally, these phases 
correspond to the first and second order expansion of the 
quasiparticle spectrum with respect to the angle formed 
by the particle momentum and the axis of spontaneous 
symmetry breaking. 

Consider a uniform gas of Fermi atoms with two hyper- 
fine states, which we assign labels 1 and 2 (these states 
equivalently can be thought of as pseudospins f and j.) 
The model Hamiltonian that describes our system is 

H = ^e p a p(J a pa - g x ^_ p , 2 a- Pt2 a p s, (2) 

p : a pp' 

where a) p(T and dp(j are the creation and annihilation oper- 
ators of a state with momentum p, pseudospin <r(= 1, 2) 
and energy e p — p 2 /2m, where m is the atom bare mass 
(here and below we set the volume V = 1). The two-body 
coupling constant is determined by the s-wave scattering 
length a < as g = AnTi 2 \a\/m. In the following we 
will work in a scheme where the particle number con- 
servation is explicitely implemented by fixing their den- 
sities pi( 2 ) = Dp""P,<T; "hp,* = Op.i(2)"p4(2) ( or equiva- 
lently the total density p and the asymmetry parameter 
a) by adjusting the chemical potentials of the hyper- 
fine states. There is a significant difference between the 
scheme above and the one where the total chemical po- 
tential is fixed and the gap is studied as a function of the 
difference in the chemical potentials of species; in the lat- 
ter case double valued solutions appear which are absent 
in the former case [l2T |. 

The mean-field solutions for the model Hamiltonian 
(|2*|) can be obtained by diagonalizing the Hamiltonian 
with the help of the familiar Bogolyubov transformations: 
bp.i = u p a p ^+v p a}_ p 2 and b Py2 — u p a p ^—v p a}_ pl , where 
u p + v p = 1. A variational minimization of the energy 
with respect to the parameter u p (or v p ) leads to the gap 
equation 

A = gjj^u p v p [l- f{Eh) f(E 2 )} , (3) 

where f(E) = [1 + exp(i?/T)] _1 is the Fermi distribution 
function, T is the temperature. The two branches of 
quasiparticle spectra are defined as 

% } = ^Je s + A 2 ±u, (4) 

where the symmetrized £s = ■k(si+£2) and anti- 
symmetrized £a — \{ e i~ e 2) spectra are written in 



terms of the normal state spectra e CT = e p — \i„ (we do 
not distinguish the masses of different hyperfine states) 
and the transformation parameters are defined as 

The occupation of the states in the superfluid phase are 
given by 

n Pjl ( 2 ) = {u 2 p f(E m ) + v 2 p [l - f(E m )]} , (6) 
with the normalization condition 

P° = ^2 n p.<?- ( 7 ) 
p 

We now turn to the description of the perturbations of 
the Fermi surfaces from the spherically symmetric form 
and the study of the stability of these perturbations. The 
two Fermi surfaces in momentum space are defined by the 
equations e a — e p — fx a = 0. When the chemical poten- 
tials [i a = p 2 Fa /2m, where pf,o are the Fermi-momenta 
of the hyperfine states, are isotropic in the momentum 
space the Fermi surfaces are spherical. Relaxing the lat- 
ter assumption we expand the quasiparticle spectrum in 
spherical harmonics e a = ei a Pi(x), where x is the co- 
sine of the angle formed by the particle momentum and a 
randomly chosen symmetry breaking axis, Pi(x) are the 
Legendre polynomials. The I — 1 terms break the trans- 
lational symmetry by shifting the Fermi surfaces without 
deforming them; these terms are ignored below. Truncat- 
ing the expansion at the second order (1 = 2), we rewrite 
the spectrum in a form equivalent to the above one 

£<r = £p - AV (l + V<tX 2 ) , (8) 

where the parameters r\ a describe the quadrupole defor- 
mation of the Fermi surfaces. It is convenient to work 
with the symmetrized and anti-symmetrized combina- 
tions Se = (771 — T)2)/2 and 5 = (iji + 772 ) / 2 , which we will 
refer to as 'relative' and 'conformal' deformations. Our 
next task is to examine the energy of the superfluid state 
at finite deformations to assess whether the deformations 
lower the energy of the system and lead to a new stable 
ground state. We shall work at fixed temperature and 
number density of the hyperfine states and will exam- 
ine the difference between the free-energies of the super- 
fluid state with deformations and the undeformed normal 
state. We shall assume that the conformal deformation 
is absent 3 = and look for a minimum of this difference 
with respect to a single parameter Se. (The situation is 
similar to the description of the LOFF phase where one 
allows for the normal state spectrum to have a finite mo- 
mentum P and seeks the minimum of the appropriate 
thcrmodynamical potential as a function of P). 

The free-energy of the superfluid phase is defined as 

Fs — ^kin + Spot — TSs, (9) 
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where the first two terms comprise the internal energy 
which is the statistical average of the Hamiltonian 
and Ss is the superfluid entropy, defined by the well- 
known combinatorial expression. In our model the sum 
of the kinetic and potential energies is 



Skin + Spot = J! £ P ("'P' 1 + "'P' 2 ) ' ( 10 ) 

P 5 

The free-energy of the undeformed normal state follows 
by setting in the above expressions A = = Se. Because 
of the contact form of the interaction the gap equation 
and the superfluid kinetic energy need a regularization. 
The regularized gap equation is 



1 = 
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dpp / dx 



l-f(E 1 )-f(E 2 
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A 2 e pJ 

(11) 

where the case 7 = 1 and A — > 00 corresponds the 
common practice of regularization which combines 
the gap equation with the T-matrix equation in the free 
space. The case 7 = and finite A corresponds to the 
cut-off regularization of the original gap equation. The 
term -Ekin is regularized by a cut-off A, which is deduced 
from the requirement that both regularization schemes 
give the same result for the gap. 

Eq. Ill 1[1 was solved numerically with the constraint Q 
for various values of the dimensionless parameter kpa at 
density p = 3.8 x 10 12 cm -3 and temperature T = 10 
nK. This density corresponds to kp rj 4.8 x 10 4 cm -1 
and Fermi-temperature Tp = 942 nK. The triplet scat- 
tering length in vacuum for 6 Li atoms in the hypcrfinc 
states |1) = \F = 3/2, m F = 3/2} and |2) = [3/2, 1/2) 
is a = — 2160as, but as already pointed out can be easily 
manipulated using Feshbach resonances. Figure ^ dis- 
plays the dependence of the pairing gap and the free- 
energy difference AF = F$ — Fn on the relative defor- 
mation for several density asymmetries and zero confor- 
mal deformation. We restrict the density asymmetry to 
positive values, i.e. assume p\> pi\ positive values of Se 
correspond to a prolate (cigar-like) deformation of the 
majority and oblate (pancake-like) deformation of the 
minority population's Fermi-spheres; for negative Se the 
reversed is true. When there are no deformations, the an- 
tisymmetric part of the quasiparticle spectrum Q, £4, 
acts in the gap equation J3J to reduce the phase space 
coherence between the quasiparticles that pair; (when 
£4 = the BCS limit is recovered with equal occupa- 
tions for both particles and perfectly matching Fermi 
surfaces). This blocking effect is responsible for the re- 
duction of the gap with increasing asymmetry and its 
disappearance above a ~ 0.07 [see Eq. JJ)]. Allowing 
for deformations introduces a modulation of £4 with the 
cosine of the polar angle x (in the frame where the z axis 
is along the symmetry breaking axis), which acts to re- 
store the phase space coherence for some values of x at 
the cost of even lesser coherence for the remainder val- 
ues. The result, seen in Figure ^ is the increase of the 




FIG. 1: (Color online) The dependence of the pairing gap 
(upper panel) and the free-energy difference (lower panel) on 
relative deformation for kpa = 0.55, at temperature T = 10 
nK, density p = 3.8 x 10 12 cm" 3 , and constant a — 0.0 (solid 
line), a = 0.02 (dashed line), a = 0.04 (dashed-dotted line), 
a = 0.05 (dotted line) and a = 0.057 (short-dashed line). 



gap for finite deformations. At extreme large asymme- 
tries the re-entrance effect sets in: the gap exists only 
for the deformed state, with lower and upper critical de- 
formations marking the pairing regions. These features 
are seen for both the positive and negative values of Se, 
but are more pronounced for Se > 0; (quite generally 
our equations are not symmetric under the interchange 
of the sign of deformation, except when a = 0). The 
free-energy difference AF mimics basically the gap func- 
tion due to the contribution from the potential energy; 
note that the critical values of deformations at which 
AF vanishes do not coincide with those for the gap due 
to the positive contribution of the kinetic energy differ- 
ence. The same calculations as above were carried out for 
larger couplings kpa — 1 and kpa = 2 with qualitatively 
similar results; the gaps found in the symmetric case are 
1.93 and 3.75 nK, respectively, the reentrance effect is ob- 
served in each case for asymmetries around 0.18 and 0.3 
and the pairing disappears above the asymmetries 0.22 
and 0.43. Figure [21 compares the quasiparticle spectra 
Ei and E% for combinations of a and Se. An important 
feature of the asymmetric (a 7^ 0, Se — 0) spectrum is 
its gapless nature, i.e. the existence of nodes for one (or 
both) branches of the spectra (c.f. with the gapped BCS 
spectrum also shown in Fig. 0). The spectra of the DFS 
phase cover the range bounded by the curves with x = 
and x = 1. We conclude that the spectrum of deformed 
superfluid states is likewise gapless for a range of the an- 
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system was in the deformed superfluid state one would 
detect a mean momentum of particles of type 1 (major- 
ity) in the direction of symmetry breaking by about 20% 




FIG. 2: (Color online) The dependence of the quasiparticle 
spectra of two hyperfine states E\ and E2 on the momentum 
for a — — Se (solid line); a = 0.05 and Se — (dashed lines); 
a = 0.05, Se = 0.1, x — (dashed-dotted) and x = 1 (short- 
dashed lines). The Fermi-momentum Rf = 4.8 10 4 cm -1 is 
indicated by the vertical line. The remaining parameters are 
as in Fig. 1. 

gles defined by the variable x. The macroscopic features 
of the atomic DFS phase, such as responses to the den- 
sity perturbations or electromagnetic probes, and their 
thermodynamic functions (heat capacity, etc) would dif- 
fer from the ordinary BCS phase due to the nodes and 
anisotropy of their spectrum, as is the case for gapless 
and/or anisotropic metallic superconductors. Figure |3| 
shows the occupation numbers in BCS, deformed and/or 
asymmetric cases; varying the cosine of the polar angle x 
covers a range of probabilities which includes the unde- 
formed asymmetric state. The bell-shaped curves show 
the angular polarization of the occupation numbers de- 
fined as 8n a — \n a (x = 1) — n a (x = 0)|. We observe 
up to 20% anisotropy in the occupation probabilities of 
particles along and orthogonal to the symmetry breaking 
axis. 

In closing, we would like to address the issue of an ex- 
perimental detection of the DFS phase. A direct way 
to detect the DFS phase is the measurement of the 
anisotropy in the momentum distribution of the trapped 
atoms. Such a measurement can be realized by the time- 
of- flight technique 0, Q • This method uses the fact that 
after releasing the trap, the atoms fly out freely and an 
image of their spatial distribution taken after some time 
of flight provides information on their momentum distri- 
bution when confined inside the trap. Assuming that the 



FIG. 3: (Color online) The dependence of the occupation 
probabilities of two hyperfine states on the momentum. The 
Fermi-momentum Rf = 4.8 10 4 cm -1 is indicated by the ver- 
tical line. The labeling of the lines is as follows: a = = Se 
(solid line), a = 0.05 and Se — (dashed lines); a = 0.05, 
Se = 0.1, x = (dashed-dotted) and x = 1 (short-dashed 
lines). The bell-shaped curves show the anisotropy - the dif- 
ference between the x = 1 and x = occupation numbers - 
for a = 0.05, Se = 0.1. The remaining parameters are as in 
Fig. 1. 



larger than that of particles of type 2 (minority) in the 
same direction. Therefore, the presence of anisotropy 
in the detected momentum distributions is an evidence 
for a deformed superfluid state being the ground state 
of the system, as deformation alone (i. e. without pair- 
ing) would not lower the energy so as to produce a de- 
formed non-superfluid ground state. The direction of 
spontaneous symmetry breaking (in /c-space and, there- 
fore, also in real space) is chosen by the system randomly 
and needs to be located in an experiment to obtain max- 
imum anisotropy. 

The Barcelona-Tubingen collaboration was supported 
in part by the programs DGICYT No. BFM2002-01868 
and HA2001-0023 [Acciones Integradas] (Spain) and the 
Deutsche Akademische Austausch Dicnst and Sonder- 
forschungsberich 381 of the Deutsche Forschungsgcmcin- 
schaft (Germany). J. M.-P. would like to acknowledge the 
participants of the Levico Workshop on Ultracold Fermi 
gases, and particularly R. Combescot and N. Nygaard, 
for valuable discussions. J. M.-P. also acknowledges the 
support from the Generalitat de Catalunya. 



5 



[1] H. T. C. Stoof, M. Houbiers, C. A. Sackett, and R. G. 

Hulet, Phys. Rev. Lett. 76, 10 (1996); 

M. Houbiers, R. Ferwerda, H.T. C. Stoof, W. I. 

Mc Alexander, C. A. Sackett, and R. G. Hulet, Phys. Rev. 

A 56, 4864 (1997). 
[2] B. DeMarco and D. S. Jin, Science 285, 1703 (1999); 

K. M. O'Hara, S. L. Hemmer, M. E. Gehm, S. R. 

Granade, and J. E. Thomas, Science Dec 13, 2179 (2002); 

C. A. Regal, C. Ticknor, J. L. Bohn, and D. S. Jin, PRL 

90, 053201 (2003). 

K. E. Strecker, G. B. Partridge, and R. G. Hulet Phys. 
Rev. Lett. 91, 080406 (2003); 

M. W. Zwierlein, C. A. Stan, C. H. Schunck, S. M. 
F. Raupach, S. Gupta, Z. Hadzibabic, and W. Ketterle 
Phys. Rev. Lett. 91, 250401 (2003); 

M. Greiner, C. A. Regal, and D. S. Jin, Nature 426, 537 
(2003). 

C. A. Regal, M. Greiner and D. S. Jin, Phys. Rev. Lett 
92, (2004) 040403. 
[3] W. C. Stwalley, Phys. Rev. Lett. 37, 1628 (1976); 

E. Tiesinga, B. J. Verhaar, and H. T. C. Stoof, Phys. 
Rev. A 47, 4114 (1993); 

Ph. Courteille, R. S. Freeland, D. J. Heinzen, F. A. van 
Abeelen and B. J. Verhaar Phys. Rev. Lett. 81, 69 (1998). 
[4] J. L. Bohn Phys. Rev. A 61, 053409 (2000); 

M. Holland, S. J. J. M. F. Kokkelmans, M. L. Chiofalo, 

and R. Walser Phys. Rev. Lett. 87, 120406 (2001); 

E. Timmermans, K. Furuya, P. W. Milonni, and A. K. 

Kerman, Phys. Lett. A 285, 228 (2001); 

M. L. Chiofalo, S. J. J. M. F. Kokkelmans, J. N. Milstein, 

and M. J. Holland Phys. Rev. Lett. 88, 090402 (2002); 



Y. Ohashi and A. Griffin, Phys. Rev. Lett 89, 130402 
(2002); 

J. Stajic, J. N. Milstein, Qijin Chen, M. L. Chiofalo, M. 
J. Holland, and K. Levin, cond-mat/0309329 
G. M. Falco and H. T. C. Stoof, Phys. Rev. Lett. 92, 
130401 (2004); 

J. Kinnunen, M. Rodriguez, and P. Torma, 
|cond-mat/0401543l 
[5] C. A. Regal, M. Greiner, and D. S. Jin, PRL 92, 040403 
(2004). 

[6] M. W. Zwierlein, C. A. Stan, C. H. Schunck, S. M. 
F. Raupach, A. J. Kerman, and W. Ketterle, PRL 92, 
120403 (2004). 

[7] J. Kinast, S. L. Hemmer, M. E. Gehm, A. Turlapov, and 

J. E. Thomas, PRL 92, 150402 (2004). 
[8] J. Mur-Petit, A. Polls, and H.-J. Schulze, Phys. Lett. A 

290, 317 (2001); 

U. Lombardo, P. Nozieres, P. Schuck, H.-J. Schulze, and 

A. Sedrakian, Phys.Rev. C 64, 064314 (2001). 
[9] H. Mther and A. Sedrakian, Phys. Rev. Lett. 88, 252503 

(2002); Phys. Rev. C 67 015802 (2003). 
[10] A. I. Larkin and Y. N. Ovchinnikov, ZhETF 47, 1136 

(1964) [Sov. Phys. JETP 20, 762 (1965)]; 

P. Fulde and R. A. Ferrell, Phys. Rev. A 135, 550 (1964). 
[11] R. Combescot, Europhys. Lett. 55, 15 (2001). 
[12] A. Sedrakian and U. Lombardo, Phys. Rev. Lett. 84, 602 

(2000); 

A. Sedrakian, T. Aim, and U. Lombardo, Phys. Rev. C 
55, 582 (1997). 



